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Abstract 

We consider an interacting particle system (r;t) t >o with values in {0, 1} Z , 
in which each vacant site becomes occupied with rate 1, while each con- 
nected component of occupied sites become vacant with rate equal to its 
size. We show that such a process admits a unique invariant distribu- 
tion, which is exponentially mixing and can be perfectly simulated. We 
also prove that for any initial condition, the avalanche process tends to 
equilibrium exponentially fast, as time increases to infinity. Finally, we 
consider a related mean-field coagulation-fragmentation model, we com- 
pute its invariant distribution, and we show numerically that it is very 
close to that of the interacting particle system. 

Key words : Stochastic interacting particle systems, Equilibrium, Coalescence, 
Fragmentation, Self organized criticality, Forest-fire model. 

MSC 2000 : 60K35. 

1 Notations and main results 

Consider an independent family TV = ((N t (i))t>o)iez of Poisson processes with 

rate 1. In the whole paper, such a family will be called an IFPP. 

Assume that on each site i £ Z, snow flocks are falling according to the process 

(Nt(i))t>o- When a flock falls on a vacant site of Z, this site becomes occupied. 

When a flock falls on an occupied site i £ Z, an avalanche starts: the whole 

connected component of occupied sites around i becomes vacant. 

We denote by (r]t(i))t>o,iez the process defined, for t > and i £ Z, by rj t {i) = 1 

(resp. rjt(i) = 0) if the site i is occupied (resp. vacant) at time t. 
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To avoid infinite rates of interaction, we will restrict our study to the case where 
the initial condition 770 lies in the following space: 

E := irj g {0,1} Z ; liminfrj(i) = liminf rj{i) = i . (1.1) 

I i — > — 00 i — > + oo j 

A state 77 belongs to E if and only if it has no infinite connected component of 
occupied sites. This condition is not really restrictive: easy considerations show 
that even if 770 € {0, 1} Z \E, r]t S E for all t > 0. This comes from the fact that 
infinite connected components of occupied sites have an infinite death rate. 

It is standard and easy to show, using for example a graphical construction, that 
for any initial condition 770 £ E, for any IFPP N, the process (j]t)t>o exists, is 
unique, and takes its values in E. It is actually a deterministic function of 7/0 
and (N(i))i£z. We call this process the 770-avalanche process, or the (770, AO- 
avalanche process when this precision is needed. See [15] for many examples of 
graphical constructions. 

Furthermore, the process (r)t)t>o is a strong Markov process, and its infinitesimal 
generator A is defined, for 77 G E and $ : E 1— ► M sufficiently regular (e.g. $(77) 
depending only on a finite number of coordinates of 77) by 

A$(r))=Y,M<v))-Hri)}, (1-2) 

iGZ 

where 0,(77) G E is defined in the following way: 

• if 77(7) = 0, then [aj(77)](i) = 1 and [ai(r))](k) = rj(k) for all k ^ i; 

• if 77(7) = 1, set l v (i) = sup{k < i;r/(k) = 0} + 1, r. q {i) = inf{fc > i;r)(k) = 
0} — 1, and put [0^(77)] (k) = for k e [Z,j («) , 7*^ («)] and [ai(r))]{k) = r](k) 
for all k £ [ln(i), r n (i)]. 

Our main result in this paper concerns the invariant distribution of the avalanche 
process. 

For A C Z and T e V(E) we denote by T A = Top^ g V({0, 1} A ) its restriction 
to A, where pa ■ E *—> {0, 1} A is the canonical projection. 

For two probability measures /a, v on a measurable space (F,f), we denote by 
\fi — v\tv = sup Ge:F |/i(G) — v{G)\ the total variation between /_t and v. 

Theorem 1.1 (a) The avalanche process admits an unique invariant distribu- 
tion n g V(E). 

(b) The exponential trend to equilibrium holds in the following sense. For if G E, 
denote by Itf the law of the ip-avalanche process at time t. There exist some 
constants C > 0, a > such that for all t > 0, all I > 0, 

sup I (nf) t _ M] - n t _ M] I < £7(1 + l)e- at . (1.3) 

(c) For I > 0, there exists an explicit (see Appendix\^j and perfect simulation 
algorithm for a ILjn- distributed random variable (?7o(*))ie[— z,/] ■ 
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(d) The invariant distribution II is exponentially mixing in the following sense: 
one may find some constants C > 0, < q < 1, such that for any k G Z, 
n G Z + , 

|n(_ 00) fc]u[fc+n,ao) _ ^-(-oo,k] ® n[fc+n,oo) | Ty < Cg™. (1.4) 

Let us comment on these results. First, the system is very stable, in the sense 
that no large clusters of occupied sites may appear. Indeed, large clusters have 
a large death rate. Clearly, the existence of invariant distributions should easily 
follow from such an argument. Of course, uniqueness of the invariant distribu- 
tion and trend to equilibrium are not suprising, but much more work is required, 
especially to give a rate of convergence. The perfect simulation algorithm we 
give is quite complicated, but gives, in some sense, an explicit expression of 
the invariant distribution II. Finally, point (d) explains that at equilibrium, for 
two sites i and j, the dependance between rj(i) and r](J) decreases exponentially 
fast with \i — j\. Such a result is also quite natural, but the proof is quite 
complicated. 

At the end of the 80's, the so-called self-organized critical (SOC) systems became 
rather popular. They are simple models supposed to enlight temporal and spa- 
tial randomness observed in a variety of natural phenomena showing long range 
correlations, like sand piles, avalanches, earthquakes, stock market crashes, fire 
forest, shape of mountains, of clouds, ... Very roughly, the key idea (present in 
Bak-Tang-Wiesenfeld [2] about sand piles) is that of systems growing toward a 
critical state and relaxing through catastrophic events (avalanches, crashes, fire, 
...) ; if the catastrophic events become more and more probable when approach- 
ing the critical state, the system spontaneously reaches an equilibrium close to 
the critical state. 

SOC systems commonly share other features as long range correlations, power 
laws for the amplitude of catastrophic events, spatial fractality of observed pat- 
terns, lack of typical scale, ... The most classical model is the so-called sand 
pile model introduced in 1987 in [2J, but a lot of variants or related models 
have been proposed and studied more or less rigorously, describing earthquakes 
(Olami-Feder-Christensen, [16]) or fire forest (Henley [13] ; Drossel-Schwabl, [8]) 
to mention a few. For surveys on the subject, see [7] or [3], for instance. 

Initially, our process was thought as a very rough simplification of a sand pile 
model. In sand piles geometric rules describe the structure of a stable sand pile. 
Sand grains fall on a given pile ; if the new pile is unstable, it is re-organized 
to become stable, through (possibly many successive) elementary steps ; such 
events are called avalanches. If the pile lives on a bounded domain, grains falling 
out of the domain disappear ; if the model is realistic, one can imagine that the 
number of grains in the pile and the shape of the pile reaches an equilibrium. 
Frequence and amplitude of avalanches at equilibrium are related to the num- 
ber of grains that disappear. In our much simpler model, a grain falling on an 
occupied site yields an avalanche involving all grains in the connected compo- 
nent (that immediately disappear). It does not pretend to be a good physical 
description of a sand pile: the purpose is more to catch what is really important 
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in SOC systems. 

This simplification is pertinent in that it can also be viewed as a particular case 
of a maybe more natural simplification of forest fire models. Roughly fire forest 
model can be described as follows: on a lattice, trees are born (sites become 
occupied) at a certain rate, say 1 ; at each tree, a fire may start at some rate, 
say A > 0: the site becomes vacant and fire propagates to neighbouring trees 
(occupied sites) at a given speed (see [5] for a precise description). Taking an 
infinite propagation speed means that the whole connected component (of sites 
occupied by trees) containing the ignited tree burns at once (one may think 
of lightning). Our model corresponds to the case A = 1, infinite propagation 
speed and a lattice equal to Z. From the point of view of SOC systems, the 
interesting phenomenon is in the asymptotic regime A — > 0. Indeed fires are 
less frequent, but when they occur, destroyed clusters may be huge. These 
models have been subject to a lot of numerical and heuristical studies (see [12] 
for references), but fewer rigorous results. Even existence and uniqueness of 
the process for a multidimensionnal lattice and given A has been proved only 
recently [HI IIP) . Limiting rescaling when A — > has been studied numerically 
[51112] but attempts to give a rigorous basis, even in dimension 1 are more recent 
[51I31[S]. Still our model had not received a complete rigorous treatment, and 
as far as we understand, even if results are not surprising they are now quite 
complete and the approach we propose may be extended. 

Consider the model in which birth flocks follow Poisson processes with rate 1, 
while killing flocks follow Poisson processes with rate A > 0. We believe that 
our result could be extended without difficulty to the case where A > 1 (so 
that the clusters are not very large). In the case where A < 1, the method we 
use probably breaks down, but the refined version of the algorithm described in 
Appendix [A] gives hints for further research in this direction 

The paper is organized as follows. In Section we show that the avalanche 
process can be coupled with (and compared to) a very simple system of inde- 
pendent particles which we call a Bernoulli process. The invariant distribution 
of this particle system is an infinite product of Bernoulli distributions. 
In Section [3J we show how to build the invariant distribution of the avalanche 
process from a stationnary Bernoulli process on an a.s. finite time intervall, 
provided some cluster (concerning essentially the Bernoulli process) is a.s. finite. 
We obtain some large-deviation type upperbounds for the width and height of 
this cluster in Section [U 

This allows us to conclude the proof in Section O the invariant distribution 
exists and can be perfectly simulated. We can estimate the decay of correlations 
in the invariant distribution of the avalanche process, using the upperbound of 
the width of the previously cited cluster. The coupling also shows, in some 
sense, the uniqueness of the invariant distribution and the trend to equilibrium. 
The rate of return to equilibrium is obtained as a corollary of the upperbound of 
the height of the cluster. In Appendix [Al we write down the perfect simulation 
algorithm for the invariant distribution derived from Sections [3J and Q] 
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We finally introduce a related coagulation-fragmentation mean-field model in 
Section [6] assuming that the correlations between the sizes of connected com- 
ponents of occupied sites are neglictable, we write down an infinite system of 
ordinary differential equations satisfied by the concentrations of clusters with 
size k, for k > 1: each pair of clusters coalesce at constant rate, while each 
cluster with size k breaks into clusters with size at rate k. The equilibrium 
state of the system of O.D.E.s can be computed almost explicitely. Numerical 
experiments show that this model is an excellent approximation of the avalanche 
process, at least from a global point of view. 



2 The coupling with a Bernoulli process 

The starting point of our results is that we may deduce a realization of the 
(possible) equilibrium II of the avalanche process from that of a much simpler 
process, which we now describe. 

Consider as before an IFPP N, and an initial state Co € {0, 1} Z . Assume that the 
snow flocks are falling on each site i according to N(i), but that the avalanche 
is restricted to the site i: if i was vacant, it becomes occupied as before, but 
if it was occupied, it becomes vacant, letting its neighbors enjoy their own life. 
Denote, for each igZ, each t > 0, by — 1 (resp. = 0) if the site i is 
occupied (resp. vacant) at time t. 

The process (Ct)t>o is obviously well-defined, unique, and explicit: for i 6 Z, 
t>0, (t{i) = Co{i) (resp. (t{i) = 1 — Co(*)) if N t (i) is even (resp. odd). In other 
words, 

Ci (i) = ![!-(- l)CoW+iv tW] , (21) 

We call it the Co-Bernoulli process (or if necessary the (Co, -/V) -Bernoulli process). 
Let us now describe its trend to equilibrium. 

Lemma 2.1 Let T — {h$o + 5^1) be the infinite product of Bernoulli laws 
with parameter 1/2. For Co G {0, 1} Z , denote by T^ the law of the ^-Bernoulli 
process at time t. Then for all I > 0, all t > 0, 



sup 

Co6{0,l} z 



(It ) H)I] -iy M] < (2l + l)e- 2t . (2.2) 

TV 



As a consequence, T is the only invariant distribution of the Bernoulli process. 

Proof Let thus Co £ {0, 1} Z be fixed and (Ct)t>o be the (Co, iV)-Bernoulli process, 
for some given IFPP N. First of all observe that for any t > 0, any i € Z, using 
the explicit expression of Ct(*) leads us to 

P[Ct(i) - 0] = l {Co(i)=0} P[iV t (i) is even] + l {Co(i)=1} P[iV t (i) is odd] 

= ^[l{CoW=o}(l + ^ 2t ) + l{CoW=i}(l-e^)]. (2.3) 
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This implies that |P[Ct(«) = 0] - \\ < er 2t /2. By the same way, \P[( t (i) 
1] - \\ < e- 2t /2, so that we get (T c t °) {i} - T {i} 



< e 21 . The result follows, 

TV 



since [—1,1] contains 21 + 1 sites and since the coordinates of ((t(i))ie[-i,i] are 
independent. □ 

Next, we explain how to reverse time in the stationnary Bernoulli process. This 
will be usefull to build from the past the invariant distribution of the avalanche 
process. 

Lemma 2.2 Let ( be a {0, l} z -valued random variable with law T. Consider 
N an IFPP, and let (Ct)t>o be the (stationnary) (Co, N)- Bernoulli process. 
Consider the process (Ct)te(-oo.o] defined by Q = C(-t)- f or a U t < 0. Then 
this is again a (stationnary) Bernoulli process, in the sense that for any T < 0, 
the process (Ct)te[T.o] * s a (0r> N T )- Bernoulli process with (r independent of 
N T , with Ct ~ r, and where the IFPP N T on [T, 0] is defined by N?(i) = 
N (-T)-(i) ~ N (-t)-(i) forte [T,0] audi G Z. 

Proof Let T = — S < be fixed. Using the explicit formula, we know that for 
all i e Z, alH € [T,0], 

Ct(i) = C(-* ) -W = i[i-(-i) CoW+Ar( - t '- (4) ] 

= I[l _ (_l)Co«+iV ( _ T) _(i)-iV t T (i)] = I[i _ (_ 1 )CT( l )+A f t T («)] j (2.4) 

since one easily checks that ( n (i) + Nr_ T \(i) and Ct(«) = C(-t)-(*) have the 
same parity. Thus, we just have to prove that (i) (t ~ I\ (ii) N T is an IFPP 
on [T, 0], (hi) Ct and 7V T are independent. 

Point (i) is obvious from the stationnarity of ((t)t>o, while point (ii) is a well- 
known fact about Poisson processes. To prove point (iii), it suffices to notice 
that for any i E Z, x E {0, 1}, 

P [Ct(0 = a H(AfW)te[T,o] ] - P [&-(0 = * HW))te[o,s]] 
= a {N s _ (l) is even} p [Co(0 = a;] + l {JVs _ (i) is odd} P [Co(*) = 1 - a:] 

Ct(0 = x] . (2.5) 
This ends the proof. □ 
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We will also need later the following monotonicity result about the Bernoulli 
process. 

Lemma 2.3 Consider N and V two independent IFPPs. Let Cd>Co G {A 1} Z - 
Consider the (Q, N) -Bernoulli process (C+ )t>o- 

T/iere exists M an IFPP such that, denoting by (C* 2 )t>o the (£q,M)- Bernoulli 
process, a.s., for all t >0, all i E Z, 



G 



(i) M t (i) = J * U {c i_ (i)=C 2_ (i)} diV a (i) + J * l {c i_ (i) ^ C 2_ (i)} dV s (i); 

(ttj 1/7* := mf{t > 0; (*) = C?{i)}, P[li >t}< e~ 2t , and (M 74+t - M 7i ) t > = 

(JV 7i+t -JV 7i ) t > ; 

(in) ifCl(i) = C t 2 (i) «ien tf+.W = tf +8 (0 M of/ s > 0; 

(w) «/CoW < Co W; then CtW < C?W a - s - f or all t > 0. 

We will say that (CtiC 2 )t>o are the {(o > Co >N,V) -coupled Bernoulli processes. 

Of course, the more natural coupling consisting in building the two Bernoulli 
processes with the same IFPP would not preserve order as time evolves. 

Proof The coupling we use here consists in choosing the same Poisson process 
N(i) for both processes when £o(*) = Co W> so t na t they will appear or die 
simultaneously, and will remain equal for all times. But if = Co W < Co W = 1 
(resp. = Co W < CoW = 1)' we use nrs ^ independent Poisson processes: 
dies using V t (i), while CtW appears following N t (i): after this first jump, they 
become equal, and we then use the same Poisson process N t (i), and they remain 
equal for all times. 

More rigorously, for i £ Z, denote by 7$ (resp. Si) the first instant of jump 
of N(i) (resp. V{i)), and put r» = Ti A S». It is immediate that r, follows an 
exponential distribution with parameter 2. Define the process M(i) by 

M t {i) = l {c i (i)=c g (i)} JV t (i) + l {cWcg(i)} [TWO + W ~ JV T4 )l{t>r 4 }] • 

(2.6) 

Then M(i) is classically a Poisson process with rate 1. We thus may define the 
(Co , M)-Bernoulli process (C 2 )t>o- 

Let us check points (i), (ii), (iii), and (iv). If CdW = Co(*)> these points are 
obvious and 7, = 0, since then M(i) — N(i) and Ct"(*) = C* 2 (*) f° r times. 
If = Co(i) < Co 2 W = 1. thcn 7* = n and = Ct (0 < C 2 (*) = 1 for t e 
[0,7,). Easy considerations show that for i > 7», Mt(i) and iVt(i) have an 
opposite parity, which implies that CtW = C 2 (*)- This shows points (ii), (iii), 
and (iv). Since 7$ = inf{t > 0, CtW = C 2 (*)}i point (i) can be written as 
M t (i) — Vt A7i + (N t — iV 7i )l{ t>7 .}, which achieves the proof. □ 

We now describe the coupling between the avalanche and Bernoulli processes. 

Proposition 2.4 Consider N and V two independent IFPPs. Let 770 £ E and 
Co £ {0, 1} Z . Assume that for all i £ 7L, r]o(i) < Co(*)- Consider the (Co,N)- 
Bernoulli process (Ct)t>0- 

There exists M an IFPP such that, denoting by (r/t)t>o the (r]o, M)- avalanche 
process, a.s., for all t > 0, all i £ Z, 

(i) Vt(i)<(t(i); 

(ii) M t {i) = f* l{ Va _(i)=t a _(i) } dN a (i) + f* t {Vs _ (i)<Ca _ m dV a (i). 

We will say that {(t,Vt)t>o * s the (Co, Vo>N, V)-coupled Bernoulli- avalanche pro- 
cess. 

Again here, building the Bernoulli and avalanche processes with the same IFPP 
would not preserve the order. 
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Proof The coupling is the following. For each i G Z, at time t > 0, we use: 

(a) the same IFPP N t (i) to make appear a flock in 77 and £ if 77 t _ (i) = ( t - (i) = 0; 

(b) the same IFPP N t (i) to make die the flock at i (in () or the whole connected 
component of flocks around i (in 77) if r)t—(i) = Ct-(i) — 1; 

(c) the IFPP AT t (i) to make die the flock at i (in £) and the independent IFPP 
Vt(i) to make appear a flock (in 77) if = r)t-(i) < Ct-(0 = 1- 

This construction guarantees that for all t > 0, all i G Z, 7j t (i) < CtW- The 
rigorous proof is similar to that of Lemma HOI □ 

This coupling is illustrated by Figure 1, and can be represented graphically in 
the following way. 

Graphical construction 2.5 (a) Initially, each site of 1 is occupied or not 
according to £0 or r/o. We draw black (resp. grey) segments to represent the 
marks of N (resp. V) above each site ofZ. 

(b) Next, we deduce the Bernoulli process 

when an occupied site encounters a black mark, it becomes vacant; 
when a vacant site encounters a black mark, it becomes occupied. 
(The Bernoulli process is not concerned with the grey marks). 

(c) Finally, we deduce the avalanche process 77: 

when an occupied site, say i, encounters a black mark, this makes become vacant 
the whole connected component of occupied sites around i; 
when a vacant site (say the site i, at time t) encounters a black mark, it becomes 
occupied (and so does it in the process Q) if and only if the Bernoulli process 
satisfies Ct-(i) = 0; 

when a vacant site (say the site i, at time t), encounters a grey mark, then it 
becomes occupied if and only ifCt-(i) = 1- 

This graphical construction is possible because rjo G E, which guarantees that 
for any T > 0, there are a.s. infinitely many sites i for which 770 (i) = Nt(i) = 
Vr(i) = 0, and since such sites cut the interactions. 

An immediate consequence of Proposition ^. 41 is the following, which we will use 
later. 

Corollary 2.6 Let IT be an invariant distribution of the avalanche process. Re- 
call that r is the invariant distribution of the Bernoulli process. Then II is 
stochastically smaller than T. This implies that for any random variable £ ~ r, 
we may find a random variable 77 ~ II such that a.s., for all i £ Z, rj(i) < C(*)- 

Proof First, Supp II C E, since the rate of death for each site is bounded 
below by 1. Consider i] = ( ~ n. Using Proposition 12.41 consider a 770- 
avalanche process and a (^-Bernoulli process such that a.s., for all t > 0, i £ Z, 
Vt{i) < Ct(0- Of course, T)t ~ II for all t, while ( t goes in law to T as t tends 
to infinity, due to Lemma [2.11 We conclude that for any 7 G {0, 1} Z , setting 
F T = {a 6 {0,f} z , V i G Z, a(i) > 7(1)}, II(F 7 ) < T(F y ). This says exactly 
that II is stochastically smaller than T. □ 



8 



Figure 1: Coupled avalanche and Bernoulli processes. 
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The Bernoulli (resp. avalanche) process is represented in black on the right (resp. in grey on 
the left) of each site. Initially, the Bernoulli (resp. avalanche) process is occupied on the sites 
1, 2, 3, 4, 5 (resp. 1, 2, 3, 4), and vacant on the sites 0, 6 (resp. 1, 5, 6). The Bernoulli process is 
easily constructed from the black marks: each time a site encounters a black mark, its state 
changes. Next, we have to build the avalanche process. The sites 2 and 4 are not affected 
by the first grey marks, since they are occupied. On the contrary, 5 becomes occupied when 
it encounters its first grey mark, since it is vacant and the Bernoulli process is occupied (at 
this time on this site). Next, the site 2 encounters a black mark, which kills him and its 
whole connected component of occupied sites, that is 1,2,3,4,5. Next, the (vacant) site 5 
encounters a black mark, but it does not become ocuupied, because the Bernoulli process is 
occupied. Next, the site 3 encounters a grey mark: since it is vacant and the Bernoulli process 
is occupied, it becomes occupied. But it is killed again by the site 4, which becomes vacant 
because it encounters a black mark. And so on... 



3 Coupling the invariant distributions 

Our aim in this section is to describe a way to build the invariant distribution 
of the avalanche process from that of the Bernoulli process, using the coupling 
introduced in Proposition ^. 41 Our method is based on the ideas of the famous 
Propp- Wilson algorithm, [17], which concerns Markov chains with finite state 
space. In the sequel, we will denote, for ( £ {0, 1} Z , 

E r .= {r,£E;Vi£Z, r,(i)<C(t)}. (3.1) 

Proposition 3.1 Let V and N be two independent IFPPs, and (o ~ T (re- 
call Lemma \2.1\) . Consider the (Co, N)- Bernoulli process (Ct)t>(b an d its time- 
reversed (Ct)te(-oo.o] built in Lemma lKM 

ForT £ (—oo,0] andip£E^ T , we denote by (Q , Vt' ip )te[T,o] ^ e (Ct, <P, N T , V T )- 
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coupled Bernoulli avalanche process with N^(i) — Nr-T)-(i) — Nt-t)-(i) and 

V t T {i) = V ( _T)-(») ~ V(_t)-(i) /or i e [T,0] and i € Z. 

Observe that a.s., due to Proposition \2.4\ for all S <T <t < 0, all if G 

r?^€% and rg« (3.2) 

Denote, for each i G Z, 6y (Tie re £ E is the state with all sites vacant) 

n = sup{T < 0;V^ G %, ^(t) - (3-3) 

Assume for a moment that a.s., for all i G Z, Tj > — oo. Notice that we have 
a.s., for all i G Z. aZZ si < S2 < 0, aZZ ipi E E^ + , <p% E E^ + , 

rtf + ™®=rtf + '» , "(i). (3.4) 

We iZitts may define r)o(i) = f?o I+s '°(*) (which does not depend on s < 0). 
Then H := C(rjo) is the unique invariant distribution of the avalanche process. 

It seems that the Bernoulli process is almost unusefull in this statement. How- 
ever, it allows us to couple all the avalanche processes (with different initial 
conditions) together. Furthermore, the behaviour of Tj will be studied through 
the Bernoulli process. For example, notice that Tj = > — oo if Co(*) = 0. 
Indeed, due to Proposition 12.41 we know that for all T < 0, all ip G E^, all 

s G [T,0], riJ> v {i) < C t (i), which implies that if Co(i) = (i.e. ( (i) = 0), then 
T%' v (i) = 0. Hence Tj = and 770 (i) = 0. When CoW = 1, it is much less clear 
that Tj > —00. 

Proof We split the proof into three parts. 

Step 1. Let us first explain (f3~2|) and (|3~4|) . First, the fact that for S < T < 
and ip G , 77 T ' ¥ ' G Eg T is straightforward from Proposition 12.41 Then the 
second equality in (|3.2p follows from the construction. Next, consider i G Z, 
si < S2 < 0, <pi G JS? and € £; . Due to the definition of Tj, we 

deduce that iJq' +Si,IPi (i) = r]Q t+Sl,0 (i). On the other hand, we get from (|3.2[) 

that ^+sa' e ^ T +s an< ^ 7 ?o i+Sl '°(*) = ^o* Ti+S2 (*), the latter being equal 
to TjQ Z+S2 '° {i) due to the definition of Tj. But using again the definition of Tj, we 
deduce that rfc +S2 - V2 {i) = ril l+S2 - {i). This shows (33). 

Step 2. Let us now show that n = C(r)o) is an invariant distribution for the 
avalanche process. To this aim, call (r/f)t>o the O-avalanche process. Consider 
also a bounded function $ : E 1— » M depending only on a finite number of 
coordinates, say $(ry) = $((r7(fc))m< n ) for some n > 0. We will show that 
Urn/r_ > _|_ 00 i?[$(r7°)] = i?[$(?7o)], which classically suffices to conclude. 

Consider now the processes coupled as in the statement. First, E[$>(r]j,)] — 
E MVo T '°)} for all T > 0. Next, on the set Q n (T) = {V \i\ < n, n > -T}, 
$>(t]q T ' ) = $(770) a.s. Since P[fi n (T)] increases to 1 as T increases to infinity 
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(because a.s., r_„ V ... V r„ > — oo), we deduce that liniT^+oo E[$(r] Q T ' )] — 
E[$>(rjo)]. This concludes the second step. 

Step 3. Consider another invariant distribution II' of the avalanche process. 
Let T > 0, and consider, using Lemma 12.61 a random variable tpx ~ II' such 
that ipr G E^ T a.s. Consider, as in Step 2, a bounded function $ : £7 i— > K 
depending only on a finite number of coordinates, say = $((77(&))|fc|<n) 

for some n > 0, and set f2 n (T) = {V |i| < n, n > -T}. Then on ft n (T), 
^{j%'' VT ) = $(?7o)- On the other hand, t%'' Pt ~ IT', since IT' is invariant. 
Using that P[f2 n (T)] increases to 1 as T tends to oo, we easily conclude that 
/ $cffl' = £[$(7? )]. Thus IT = n. □ 



4 The contour process 

Our aim in this section is to define and study a process which will allow us 
to estimate Ti, for i £ Z, and to bound the number of sites involved in the 
construction of ?7o(i)j m order to estimate the decay of correlations. 
The first idea is the following: consider the occupied zone in Z x [0, 00) of 
the Bernoulli process. Clearly, if this occupied zone has no infinite connected 
components, then r, is finite a.s. for all i 6 Z. Indeed, each site i would then 
a.s. be encompassed by a vacant zone of the Bernoulli process, which implies 
that the avalanche process is also vacant, and cuts the interaction in some 
sense, which would allow us to build 770 (i) from the stationnary process, using 
the graphical construction ^. 51 

But such a consideration would probably lead to a fat tail estimate of the distri- 
bution of because we are in a critical case (the proportion of space occupied 
by the Bernoulli process is 1/2). A way to overcome this difficulty is to make 
use of the grey marks (recall Figure 1), which also give us some information 
about r/o (i). 

Let us now define the left and right contour processes, keeping in mind the 
coupling between stationnary measures built in Proposition 13. II 

Definition 4.1 Let (0 6 E. and N,V be two independent IFPPs. We con- 
sider the (£0, N)- Bernoulli process (Ct)t>o o-nd we introduce the filtration Q t = 
o-{( (i),N s (i),V s (i); sG [0,t],ieZ}. 

Fori G Z, we define the (£o, N,V) -right contour process (i?J)t>o around i, with 
values in Z + ~ U {00} (see Figure 2 for an illustration) by 

n>0 

where: 

Initially, T$ = 0, i% = inf{fc > i, ( (k) = 0} - |. For n > 0, 
T* +1 = infO > Tl, ANti^ + ^ + AN^-^ + AVti^ -i) > 0}. (4.2) 
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Then 

(a) if ANjn ^Bifi + i) > 0, then R^ = inf{fc > i?^ , Cr +1 (k) = 0} - \; 

(b) if AAT T i ^{Ripi - |) > 0, then R^ = sup{fc < R l T , , Ct; +1 0) = 1} + i- 
f C J ifAVri^R^ - |) > 0, ifcera 

r»; if£Zs R ii - 1) = x ' ^t* +1 = 4«» 

(iij if( T ^ +i (R l Tl - |) = 0, t/ien i = sup{fc < J^i! Cr; +1 ( fc ) = 1} + |- 
TTie Ze/i contour process (L l t ) t >o around i is defined symmetrically. 

Remark that the sequence (T^) n >o contains all the instants of jumps of (R\) t >Q, 
but it contains also fictitious jumps (case (c)-(i)). We explain how to build 
graphically these contour processes, as illustrated by Figure 2. 

Graphical construction 4.2 Draw above each site i 6 Z the marks of N in 
black and those of V in grey. Draw in black the Bernoulli process corresponding 
to a given initial data Q$. 

A time 0, the right contour process R® lies on the left of the first vacant site of 
Co on the right of (e.g., if Co(0) = Co(l) = 1 and Co (2) = 0, then i?g = 1.5). 
Next the dynamics of R are the following: 

(a) each time it encounters a black mark on its right, it jumps to the left of the 
first vacant site on its right; 

(b) when it encounters a black mark on its left, it jumps to the right of the first 
occupied site on its left; 

(c) when it encounters a grey mark on its left (say that R®_ = i + 0.5), and if 
Ct-(i — 1) = 0, then it jumps to the right of the first occupied site on the left of 
i - 1. 

The process (L°) t >o follows the same dynamics, permuting the roles of left and 
right. 

We will see in the next section that it is possible to build rjo ~ IT (where II 
is the invariant distribution of the avalanche process) in such a way that r]o{i) 
depends only on Co, N and V in the box delimited by L\ and R\ until they first 
meet (if they do). The main reason for this is the following property, which says 
that in some sense, the contour processes encompass a given site i by a vacant 
zone of the Bernoulli process. 

Lemma 4.3 We adopt the notations of Definition \4.1\ A.s., for all t > 0, all 

iez, Ct(i4 + !) = Ct(£i-!) = o. 

Proof It is clear from the construction. □ 

Remark here that case (c)-(i) in Definition ^. H is considered to have this Lemma. 
Indeed, if we want the right contour process to have only vacant sites (of the 
Bernoulli process) on its right, we can use grey marks to jump to the left only 
when there is at least one vacant site on its strict left. 
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Figure 2: The contour processes (Fq)t>o and (L®)t>o around 0. 




The process R® (resp. L®) is represented in plain (resp. dashed) line. First, R9, = 0.5, since 
the first vacant site (of the Bernoulli process at time 0) on the right of is 1. By the same way, 
Lq = —1.5. Next, R° encounters a grey mark on its left, but since at this time (t(— 1) = 1, it 
does not jump. Then R° encounters a black mark on its right, so that it jumps to 2.5, i.e. the 
left of 3, which is (at this time) the first vacant site on its right. Next, it encounters a black 
mark on its left and jumps to 1.5, which is the right of 1, i.e. the first vacant site on its left, 
and so on... As we see on the picture, when it encouters its fourth grey mark on its left, we 
have R®_ = 3.5, and since £t-(2) = 0, it jumps to 0.5, which is the right of 0, i.e. the first 
occupied site on the left of 2. 

To study the decay of correlations, we have to estimate the width of the box, 
while to study the rate of trend to equilibrium, we have to estimate its height. 
The following estimates, central in our proof, will provide some bounds on these 
quantities. 

Proposition 4.4 Let Co ~ I\ let N, V be two independent IFPPs. Consider 
the right and left (Co, N,V)- contour processes (R^) t >o and (L®) t >o around 0. 
Consider the stopping time (for the filtration (Gt)t>o) 

p° = M{t > 0; R° t < L° t } (4.3) 

and the random variable 

iC-supi? t °. (4-4) 

*>o 

(a) There exists (3 > such that E\e$ p ] < oo. 

(b) There exists 7 > such that E[e^ R °°] < 00. 

The remainder of this section is devoted to the proof of these estimates. They 
seem quite natural, since the process (i?°) t >o is a sort of random walk with 
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negative jump size expectation: two types of events allow (R^) t >o to jump to 
the left, while only one type allows him to jump to the right. Furthermore, 
some symmetry seems to hold between the jumps to the right (due to TV) and 
those to left due to N. Thus, the result seems almost obvious, and intuitively 
very clear. However, we have not found a simple proof. Of course, the main 
difficulty is that (R®)t>o is not a continuous-time random walk: its sizes of 
jumps are not independent. Thus quite a precise study has to be done. Our 
strategy consists in bounding from above )t>o by a continuous-time random 
walk with negative jump size expectation. We first describe some immediate 
properties of the contour processes. 

Lemma 4.5 We adopt the notations of Definition ^. 1\ Let i e Z be fixed. 

(a) If Co ~ T, then the processes (i — L\) t >Q, {R\ — i)t>o, and (i?°) t >o have the 
same law (but are far from being independent). 

(b) If Rq < oo, then R\ < oo for all t > a.s. 

(c) For j < i, R? t < R\ for all t > a.s. Furthermore, if R\ = R\ for some t, 
then a.s., R 3 t+S = Rl+ S f or all s > 0. 

(d) The counting processes 

Z t' 1 = X ^{t>Tk-'\ X 1 {t>T>}l{A7V T i i (^ i +^)>o} = N t{R\- + tj) 

n>l n>l " 

Z t' % = X a {*>T^ ! } := X 1 {t>^} 1 {AAf n (i?^ -^)>0} = N t(Rt- - ^) 
n>l n>l " ™ _1 

Z t* = X 1 {t>T*>} ■= X ht>T^ { AV Tk (RU -i)>0} = Vt(Rl - \) (4.5) 
n>l n>l " 

are three independent Poisson processes with rate 1. They are (Qt)t>Q-adapted, 
and independent of Co ■ 

Remark here that Z l ' % counts the jumps to the right of R z , while Z 2 ' 1, counts 
its jumps to the left due to N (black marks on Figure 2) and Z 3 ' 1 counts its 
possible jumps to the left due to V (grey marks on Figure 2). 

Proof Point (a) is obvious by symmetry and invariance by translation. Point 
(b) follows from the fact that the Bernoulli a.s. process belongs to E for all 
t > 0, even if it does not at time 0. Point (c) is clear from the construction. 
Point (e) follows from classical properties on Poisson processes. □ 



We carry on with a natural monotonicity property. 

Lemma 4.6 We consider three independent IFPPs N , V and W . Let also 
Co, Co £ {0j 1} Z satisfy, for all i 6 Z, CoM — CoW- Then we build, recalling 
Lemma \2.3l the (Co,Co; N,W)- coupled Bernoulli process (Q, Ct)t>o- As stated 
in Lemma \2.3l (C})t>o * s the (Co, N) -Bernoulli process, while {Ct)t>o is the 
(Co, M)- Bernoulli process for some IFPP M. We denote by A ) t > (resp. 
(R°' ) t >o) the (Co 1 ,^,^) (resp. (Q,M,V)) right contour process around 0. 
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We will say that (Rf' , Rf ,2 )t>o ar ^ the (Co, Co , N, W, V)-coupled right contour 
processes around 0. We have a.s., for all t > 0, R^ 1 < i?j' 2 . 

Proof The proof is obvious from the definition of the contour process, since we 
know from Lemma |2~31 that a.s., for all t > 0, alH 6 Z, Ct (*) < (*) ■ ^ 

We consider the following initial condition. 

Notation 4.7 We say that a {0, l} z -valued random variable Co has the disti- 
bution S is if (,a(i) = 1 /ori < 0, Co(l) — 0) an d if (Co(*))i>2 *s a family of i.i.d. 
Bernoulli random variables with parameter 1/2. 

Let us now explain our strategy to bound the right contour process by a random 
walk: 

we will first upperbound the initial configuration Co of the Bernoulli process by 
a (possibly shifted) realization Co of E, we thus upperbound our contour process 
by the corresponding contour process R , 

then we will wait for the first instant ' of jump to the right of i?° to the 
right; this yields a total jump which we will call Y\ :— Rfi.o — Rq, and whose 
expectation will be shown to be negative; 

we will also observe that at T^' , we may again bound the configuration of the 
Bernoulli process by a realization Co of 5 (shifted around Rfi.o) independent of 

Yi; 

this last renewal argument allows us to build, recursively, a random walk with 
negative mean jump size, bounding from above our contour process. 

Lemma 4.8 Let Co ~ E, and consider two independent IFPPs N, V . Consider 
the (Co, N ,V)-right contour process (R^) t >Q around 0, observe that R® = 1/2, 
and denote by T^' := inf{< > 0, Ai?° > 0} the first instant where it jumps 
to the right. We also consider (Ct)t>o the (Coj N)- Bernoulli process. We set 
Yi = R° t ,, - 1/2. 

(i) ThenE\Yi] < 0. 

(ii) For all e £ (0,ln2), E[e eYl ]< oo. 

Furthermore, then there exists Co ~ E such that 

(Hi) a.s., Cfi,o(M 1)0 +i - 1/2) < for all i £ Z, 
1 1 

(iv) Co is independent ofTCfi.o, where TL t = o{R? s ,s < t}. 

Proof To simplify the notation, we omit the superscript (which says that 
we are dealing with the contour process around 0) in this proof. We consider 
the three independent Poisson processes with rate 1 (see Lemma 1431 - (d)) Z\ = 
N t (Rt- + 1/2), Z? = N t (R t - - 1/2) and Zf = V t {R t - - 1/2), and we denote 
by (3j)»>i) (Ti)i>i, (Tj 3 )i>i, respectively, their successive instants of jumps. 
We also denote by Aj = Ui>i{f/}, for j = 1, 2, 3. We set Z t = Z\ + Zf + Zf, 
which is a Poisson process with rate 3, we denote by (Ti)j>i its successive 
instants of jumps, and we set A — Ui>i{Xi} = A\ U A^ U A3. Finally, we also 
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set for convenience T = Tq = Tp 
Y 1 = R fl - 1/2. 



q — Tq = 0. Recall that we want to study 



Step 1. For n > 1, the event fi„ = {Ti (/ A 1; ...,T„_i ^ A 1 ,f n e A n } occurs 
with probability p n := 2 3 „ , since Ai,A2,A% arc the sets of jumps of three 
independent Poisson processes with same rate. Notice also that on fi„, !Zf = T„, 
and we may write Y± = — ^7=1 ^-m where X\, ...,X n are the successive 

sizes of the (possibly fictitious) jumps of R (at the instants T\ < ... < T n ), with 
X 1 > 0, X n > 0. We obtain 



E\Yi] = ]T e[{-{x 1 + ... + X„_!) + X n }l Q J. 



(4.6) 



i>i 



Step 2. Let us now bound from below C^„ := £7[JQ1qJ, for 1 < i < n — 1. 
We denote by the number of vacant sites of the Bernoulli process on the 
strict left of Rf._ — 1/2 at time , that is 



Zi := Rf - 3/2 - sup{j < R f _ - 3/2, Cf _(i) = !}• 



(4.7) 



Then, due to the definition of i?, we know that 

(a) on Tf := {f, e A 2 }, X, = = 1 + Z i; 

(b) on if := {T G A 3 }, X; = -A^ = (1 + Zj)l {Zl > 1} . 

Observe that P[T 2 |r2„] = P[if|ft„] = 1/2, and that F 2 ,F 3 are independent of 
(Z^ Tj) conditionally to 0„. These are standard properties of Poisson processes. 
Hence, 

Ci, n = ^[(1 + Zi)ln n + (l + Zj)l { z 4 >i}ln B ] 
= S[(l + Z)l n J-ip[Z = 0,n„] 

= p[n n ] + 5^p[z i >A ! ,fi n ]-ip[z i = o,n n ]. (4.8) 

fe>l 



Let now fc > 1 be fixed. We have 
P[Z > k,Q n ] =p(( fi _(R fi _ - 

= E 



|) = o,..,Cf i _(^fi_-^-*) = o,n„) 



n p (cf < -(%--^-0 = o «„,t,t_^ i n 



(4.9) 



Indeed, recalling that on 0„, R has had only jumps to the left before T, we easily 



are 



deduce that on fi n , the values of the Bernoulli process at sites j < Rf ._ 
mutually independent conditionnally to Tj,Tj_i. 

Let us set p s := (1 - e - 2s )/2 = P[N S e 2N + 1] for s > (for (AT t ) t > is a 
standard Poisson process with rate 1). 
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Now for I > 2, the site Rf._ — \ — I was occupied at time 0, and its evolution 



is obviously independent of (-Rt) te m f 

1 



so that 



P[( f (R fi _---l) = 



(4.10) 



Next, the same argument holds for I = 1 on the set > 0}, which indicates 

that the previous jump to the left was not fictitious: we have 



P[{ fi _(R fi _--)=0 



n n ,T i ,T i - 1 ,X i - 1 >0 )=Pf i . (4.11) 



But on the event = 0}, we know that Cf. 1 (Rf._ — §) = 1- Hence we get 



(4.12) 



PCf i -(%--^) = o 



Qn,Ti,Ti^i, Xi-i —01 — Pfi-fi_i 



Noting that p^. > Vf -f i> * na * = 0} C {Tj_i E ^3}, we deduce that 



(4.13) 

Gathering the estimates obtained for i > 2 and 2 = 1, we obtain, for k > 1, 
P[Z< > fcA] > E [tny^T 1 (l { f 4 _ I gA a }Pf 1 + 1 {T i _ ie ^ 3 }Pf 4 -f i _ 1 )] ( 4 - 14 ) 

Using finally classical properties of Poisson processes, we see that Tj_i,Tj are 
independent of fi n , {Ti_i G A 2 }, {T;_i G A 3 } and that P[fi„ n {T;_i G A 2 }] = 
P[Q„ n {fi_! G A 3 }] = P[Q n ]/2 = 2™- 1 /2.3", so that 



p\z t > k, n n ] > 



2.3 r 



-E 



(pft +Pf < -f i _ 1 ) P 



Ti 



(4.15) 



Next, 



P[Z t = 0, n„] = P[O n ] - P[£< > 1, O n ] < 
Thus, recalling (|4.8[) . for any 1 < i < n — 1, 



~>n — 1 



2.3" 



Pfi Pfi—Ti—i 



(4.16) 



2 n— 1 

C,;.™ > ———E 



> 



3' 

2 n-l 
^3™ 
2 n-l 
2.3™ 



-E 



PTj Pfi—fi-i ^ Pli Pfi-fi-t 

2 - 2p f . I 

Pt, Pf i -f l _ 1 , 1 
1 - Pf 1 - Pf. 



(4.17) 
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where the last equality stands for a definition. 

Step 3. We now upperbound C n , n '■= E[Xn^-Q n ]- We denote by Z n the number 
of occupied sites on the strict right of Rf n _, that is 

Z n = inf{j > R fn _ + 3/2, (f n _(j) = 0} - R fn _ - 3/2. (4.18) 

By construction, we have X n = 1 + Z n on fl n . For k > 1, we set Jk '■= 
Rf n _ + k + \ and ^ :— Cf„-(^fc)- By construction, we have X n = 1 + Z n on 
f2„, so that for k > 1 

i^n > fe, On] =P[Zl = 1, 6 - 1, fin] • (4-19) 

We now introduce the tr-field generated by the path of (Rt)te[o f containing 
also the fictitious jumps, that is, for v defined by T v — T\ {y = n on fi n ), 

W := a (i/, fi, ...,T„, A%, Afi^J . (4.20) 

Observe that obvisouly, {Jk)k>i and Q n are 7i-measurable. 
We will show that conditionnally to H, the sequence (Cfc)fc>i is a family of in- 
dependent random variables on fl n , and that for each k > 1, -P[£fc = 1|W, fi n ] < 
1/2. Since f2„ belongs to H, for all n > 1, we will deduce that 

1 1 2™ _1 

P[Z„>fc,Q„]<_P[Q„] = ___, (4.21) 

so that (since X n = 1 + Z n on 

c„,„ - ]T p[i + z n > k, n n ] = p [Zn > k, n n ] < — . (4.22) 

fc>l fe>0 

Let us thus check the announced properties of the sequence (£,k)k>i- For each 
k > 1, let Tfc = if Jfe > 2, and let be the unique instant t E [0, T-j 1 ) such that 
Jfe — 1 £ {-Rt, if Jfe < 1. We refer to Figure 3 for an illustration. Roughly, 
Tfc is the last instant before T/ where we get some information (from 7i) about 
the site Jk- Of course, (rfc)fc>i is 7i-measurable. 

We will show that the family of random variables Cr fc (Jfe) is mutually indepen- 
dent on 0„ conditionnally to H, and that for each k > 1, P[Cr fe (Jfe) = fi„] < 
1/2, This will imply the announced properties, for two reasons: 

(i) conditionnally to TL and Ct n , the evolution of the Bernoulli process at two 
different sites Jk (on [t^,^ 1 )) and Ji (on [77, T-j 1 )) arc independent, since they 
concern independent Poisson processes, 

(ii) for eachfc > 1, P[Cfi_(Jfc) - l\H,Sl n ] = P[(r k (J k ) = 1|W, 0„](l-p f i_ r J + 
P[C Tfe (Jfe) = 0\H,n n ]p f i_ Tk < 1/2. Indeed, recall that p s = (1 - e - 2s )/2 < 1/2 
stands for the probability that a standard Poisson process at time s is odd, and 
that for a, b G [0, 1/2], o(l - 6) + (1 - a)b < 1/2. 
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Figure 3: Illustration of Step 3 (and Step 6). 
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-8 -7 -6 -5=J, -3=1, -l=Ij L=J, 3=J, 

-4=X, ~ -2=| 0=J 6 2=J 8 

With this realization, we have Rfi_ = —6.5 and G = {2,4,6}. We remark that the only 
site which is occupied by the Bernoulli process when crossed by R is the site J4 = —2: it is 
crossed through a grey mark. 

Consider now the random set, measurable with respect to 7i, 

G := {k > 1; r k+1 > r k+2 } = {k>l; J k = i? Tfc - - 1/2}. (4.23) 
We notice that { J k , k S G} C [Rfi , 0]. 

We observe (see Figure 3) that conditionnally to H, f2„, for all k > 1, we have, 
( Tk (Jk) = l { J k < 1} l {keG} l {Tk eA 3 } + l {Jfc >2}Co(Jfc)- (4.24) 

Indeed, if J k > 2, the formula is obvious because then r k = 0. If J k < 1, this 
comes from the fact that the only way for R to jump to the left through an 
occupied site is that the jump follows from a grey mark (i.e. r k € A3) and that 
the concerned site is just on the left of R Tk - (i.e. k € G). 

Conditionnally to H, fl ni the sequence of events {r k £ As} ke G is independent, 
this assertion makes sense since G is itself H-measurable. Indeed, we know 
that conditionnally to Ti., for k £ G, {r k E A 3 } depends only on the Poisson 
processes N t (J k ),V t (J k ) (and on Q Tk -{Jk — 3/2) which is H-measurable). The 
conditionnal independence (with respect to H, f2 n ) of the family {r k G A^\ ke c 
follows then from the fact that for k\ < k 2 in G, J kl > J k2 (see Figure 3), and 
from the independence of the Poisson processes (N t (i), V t {i)) and (N t (j), V t (j)) 
for i^j. 

Recall now (|4.24|) . Using the conditionnal independence (with respect to H, 0„) 
of the family {r k £ A^} ke Gi the fact that the sequence (J k ) k >i is 7i-measurable, 
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and the fact that the family (Co(«))z>2 is mutually independent and independent 
of TL 1 f2„ (these are i.i.d. Bernoulli random variables with parameter 1/2), we 
obtain that conditionnally to 7i, O n , the sequence {C, Tk {Jk))k>\ is mutually 
independent. 

We finally conclude by noting that for any k > 1, 



P 



(r k (Jk) = l 



H, O n 



= i { j k < 1} t {keG} p[T k g A 3 \n,n n ] 



-Mj h >2} P 



Co(Jk) = 1 



< 1/2. 



(4.25) 



The last inequality comes from the fact that P (o(Jk) = l\H.,Q n = 1/2 if 

Jk > 2 as was previously noticed, while for k G G, P [r k G A 3 \H,h n ] < 1/2. 
Indeed, having a look at Figure 3, we realize that for k G G, {r^ G A3} C 
{( Tk -(Jk — 3/2) = 0}, and that due to classical properties of Poisson processes, 



P 



n e A 3 \n, n n , {( Tk -(j k - 3/2) = 0} = 1/2. 



(4.26) 



Step 4. Gathering (gU), (|4~TTj) and g^U), we get 



< 



I '171 ± 

E^-E^a+^H-E 



n>l 



1=1 



-(1 + Bi) 



o i — 1 

< i-eV^ 

i>l d 



(4.27) 



To conclude that < 0, we thus have to prove that / = J2i>i > 1- 



But, we may write I = Ii + I 2 + h/2 + I4/2, with 



oi-l 

h ■■= y —e 



1 ~Pfi 



02 — 1 

y — e 



(>1 



i>i 



i>i 



1 -Pfi 



PTi-fi-! 



(4.28) 



Since — is exponentially distributed with parameter 3 (for all i > 1), 



E 

i>l 



ds3e 



-3s 



-2s 



Next, since Ti follows a T(i, 3)-distribution, 



*3 = E 



ds- 



-2s 



3 /l 1 

2 \ 3 ~ 5 



dse 



-2 s 



(4.29) 



1 

= 3' 
(4.30) 
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and (using the substitution u = e s ) 

■yi — l roo ni 



* = E 

»1 







-2s 



dse 



-2s 



, 1 — It 7T 

ait ^ = 2arctanl — 1 = 1. 

1 + u 2 2 



(4.31) 



Finally, using the independence between Ti_i and Tj — Ti_i, we get 



3 



+E 

i>2 

I ^ 

3 7 



1 -Pii 

2« — 1 ^-oo 



c?s3e 



oi — 1 



II 



-3s 



(i-2)! 



dt3e' 



-3t 



1 



-2/ 



1 



-2s 



1 + e- 

1 ».2/i „,2\ 



ds 



du 



o 



w 2 (l - u 2 ) 
1 + u 2 



5 7T 

3 ~ 2 



, ^ 2 
2, 4"3 



1 

" 3' 



dt2e~ s e 



1-v 2 
1 + it 2 ?; 2 



s„-3t 



-2s— It 



-21 



(4.32) 



We finally get that 7=f — l + | + | + ^ = f — |>1. Thus #[Y X ] < 0. 

Step 5. We still have to prove the exponential moment estimate. Using the 
same notation as previously, we will just use that for each n > 1, Y\ < X n 
I + Z„ on fi„. Recalling (|4.21l) and that X)n>i = 1j we g et ; for an y fc > 1 



P[Yi > k] < E ^ > fc - 1, n»] < 



n>l 



(4.33) 



We classically conclude that for e £ (0,ln2), £J[e eYl ] < oo. 

Step 6. We finally have to build Q. First note that obviously, Cf 1 (* + Rf 1 



i/2) < i = CoW for * < o> while Cf^ 1 + ^fi - V 2 ) 







Co(l) due to 



Lemma l4~3l Hence we just have to build Co (*) f° r « > 2. We write for simplicity 
Ki = i + R f i -1/2. 

Using the same arguments as in Step 3, one may check that conditionally to 
TLfi, the family (Cfi (-?Q))i>2 is independent, and that for all i > 2, 



< 1/2. 



(4.34) 



We consider a family (f/j)j>2 of i.i.d. random variables uniformly distributed 
on [0, 1] (independent of everything else) and we set, for each i > 2, 



8(0 



1 



{Ofi(^)=i} 



1 {4 1 (K i )=0} 1 {^<<i)}' 



(4.35) 
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where, due to (|4.34|) . 



e(i) 



1 -2P[C f i(P.) = 



2P[£ f i(tf i ) = 0|W f , 
We next observe that for all i > 2, 



e [0,1]. 



(4.36) 



P 



= E 



P[%(i) = l\Hf a] 



P 



p[C f i(^) = i|w f i] + £ (i)P[C f x(^) = o|?%] 



(4.37) 



due to our choice for e(i) and to the independence of Ui of everything else. We 
deduce that for each i > 2, Q(i) is a Bernoulli random variable with parameter 
1/2, and that it is independent of Hfi. This and the conditionnal (to Hfi) 

independence of the family (Co(*))»>2 clearly imply that finally, (Co(*))i>2 is an 
i.i.d. sequence of Bernoulli random variables with parameter 1/2, independent 
of Hfi. Finally, it is clear from ([435]) that for all i > 2, > Cfi (P*)- This 
concludes the proof. □ 



The following lemma shows a way to bound from above the right contour process 
started with the initial condition Co ~ T by a continuous-time random walk. 

Lemma 4.9 Let Co ~ T, let N, V be two independent IFPPs, and consider the 
(Coi N,V) -right contour process (R®)t>o around 0. Then for k > 0, P[Rq = 
k- 1/2] = (l/2) fe+1 . Furthermore, we may find a Poisson process (Z t )t>o with 
rate 1, a family of i.i.d. random variables (5^)i>i distributed as Y\ (see Lemma 
\4-8\ ) in such a way that i?§ and ((Zt)t>o, (5^)j>i) are independent, while a.s., 
for all t>0, 

z t 

R°<R% + J2 Y i' ( 4 - 38 ) 



Proof We omit as in the proof of Lemma [4.81 the superscript 0. We consider 
Co ~ r to be fixed, and a (o-right contour process (Rt)t>o around 0. First, it is 
obvious that for k > 0, 

P[P = k - 1/2] = P[£ o (0) = 1, C (k - 1) = 0, Co(fc) = 1] = (l/2) fc+1 , (4.39) 



since Co ~ T. Next, let us explain (I4.38p . The main ideas are the following: we 
first bound our contour process by a contour process P 1 which starts from a 
(shifted) S-distributed initial data. When this process first jumps to the right, 
at some instant t\, we bound the Bernoulli process at this time by a shifted 
S-distributed data, independent of {R\ )te[o,Tx] ■ Thus we make start again a 
contour process P 2 from this S-distributcd initial data in such a way that it 
dominates (with a shift) P 1 , and thus P. And so on... The advantage of this 
method is that the increments (between two renewal times) are independent. 
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We define Co by Co(») = 1 if i < 0, CoW = 1> and Co W = Co(« + i?o - 1/2) for 
i > 2. We observe that Co(« + Ro~ 1/2) < Co W for a11 i e z > 

using Lemma 14.31 
We also notice that Co is independent of R , and is S-distributed. 
We thus may find, using Lemma a contour process (Rl) t >o around 0, start- 
ing from Co j independent of i?o, such that for all times, R t — Ro < Rj — l/2 (recall 
that R t - Ro starts from 1/2, and that for all i 6 Z, Co(« +Rq-1/2) < Co(«))- 
On the other hand, we consider t\ = inf{i > 0; AR} > 0} the first instant 
where -R 1 jumps to the right, so that R\ < R^ for all t G [0, ti]. Hence setting 
Y 1 := R^ - 1/2, we finally obtain that a.s., for all t E [0, r x ], R t <Rq + Y 1 . We 
also observe that Ro and (ri, Yi) are independent. Finally, t\ is exponentially 
distributed with parameter 1, due to Lemma 14.51 -fd), and Y\ is distributed as 
in Lemma 14.81 by construction. 

Due to Lemma 14.81 we may find Co ~ 2, independent of Rq and Yi, such that 
a.s., for all % 6 Z, Ct x (« + ^ ~ I/ 2 ) < Co where (Ct)t>o is the Bernoulli 
process starting from Co associated with the contour process (R\ )t>o- 
Using Lemma 14.61 we thus may build a contour process (R^)t>o around 0, 
independent of Ro and Yi, such that for all times R$. , t — R\ < R% — 1/2. As 
a consequence, we observe that for all t > 0, 

R n+ t ~Ro< R\+t - 1/2 < - 1/2) + (R 2 t - 1/2) <Y 1 + (R 2 t - 1/2). (4.40) 

Denote by t-i = ini{t > 0; AR 2 > 0} the first instant where R 2 jumps to the 
right, so that R 2 < R 2 T2 for all t £ [0, r 2 ]. Hence setting Y 2 := R 2 T2 - 1/2, we 
finally obtain that a.s., for all t e [ti, ti+t 2 ], R t < R +Yi+Y 2 . We also observe 
that {t2,Y 2 ) is independent of R and (t±,Yi). Finally, r 2 is exponentially 
distributed with parameter 1, due to Lemma 14.51 -fd). and Y 2 is distributed as 
Y\ by construction. 

Iterating the procedure, we find an i.i.d. family (Tfe,Y/ £ )fe> 1 of random vari- 
ables, independent of Ro, such that Y\ is distributed as in Lemma |4~51 and t\ is 
exponentially distributed with parameter 1, such that for all t > 0, 

Rt~ Ro <J2 Y ^{t>r l+ -+r k }- (4-41) 
k>l 

This ends the proof. □ 

We finally conclude the proof of the main estimates of this section. 

Proof of Proposition 14.41 We omit the superscript for simplicity. The 
proof is based on the use of Lemmas 14.91 and 14.81 We thus write, according to 
63H|, R t <R + S Zt , with S n =Yt + ... + Y n . 

First of all, we deduce from Lemma [4.81 (iVfii') that there exists 7 6 (0,ln2) 
such that q := E[e~< Yl } < 1. 

Since 7 G (0,ln2), we also deduce from Lemma l4~9l that ^[e 7 ^ ] < 00. 

Next, we recall that since (Z t )t>o is a Poisson process with rate 1, for all t > 0, 

P[Z t < t/2] < e- 5 \ where S := (1 - ln2)/2 > (any 8 > would work as well). 
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We have Roc — sup t>0 R t < Ro + sup n>1 S n . This implies that e 7 ^ 00 < 
e lR ° X)«>i e 75 ". Thus, since q £ (0, 1), and since S n and Rq are independent, 

E[e~< R ~] < E[e^ Ro ]J2 E i elSn ] < C Yl E i ejY T = C£)g n < oo. (4.42) 



n>l 



n>l 



Next we want to upperbound p. First, p < inf{i > 0, R t < and L t > 0}, so 
that by symmetry, for any t > 0, 

f[p > *] < > or L t < 0] < 2P[R t > 0]. (4.43) 

But, since R t < Ro + Sz t , 

P[p >t]< 2P[R t > 0] < 2P[Z t < t/2] + 2P[R + sup S n > 0] 

n>t/2 

< 2e~ 5t + 2E \e^ Ro+sup ™>t/i < 2e~ 5t + 2E [e lR °] ^ E [e 7S "] 



n>t/2 



< 2e' 5t + C 1 n < 2e ~ St + C 1 t/2 < Ae ~ 

n>t/2 



(4.44) 



for some constants A > 0, a > 0. We classically conclude that for any [3 G (0, a), 
£[e^] < oo. □ 



5 Proof of Theorem 11.11 

Our aim in this section is to conclude the proof of our main result. 

Proof of Theorem 11.11 We divide the proof into several steps. Let us recall 
briefly the notation of Proposition 13 . II and Definition 14. 11 we consider two inde- 
pendent IFPPs N,V, and (b ~ T. Let (Ct)t>o be the (Co, iV)-Bernoulli process 
(C*)*>o, an d let (Ct)te(-oo,o] be its time-reversed built in Lemma l2~2l 
For T G (—oo,0] and ip G E^ T (recall (I3.1|) L we denote by (Ct, ijI' V )t£[T,o] the 
(Ct, y>, -^ T , ^ T )-coupled Bernoulli avalanche process with Nf(i) = N/_T)-(i) — 
N ( _ t) _(i) and V t T (i) = V(_ T )-(i) - V(-t)-M for i 6 [T, 0] and t G Z. Recall 
that 

n = sup{T < 0;V^ G E (T , ,£>(*) = (5-1) 

We then may put rjo(i) := f?o" ,+s '°(i) (for some s < 0, recall Proposition 13. ip 
provided Tj > — oo. Recall also that by definition of t;, we have, for all T < 0, 
all G E{ T) f^ ,,p {i) — rjo(i) on the event {T < r,}. 

Next, we consider the (Co, N, V)-left and right contour processes (L\) t >a and 
(R\)t>o around i, for each i 6Z, and we adopt the notation 

R\ = sup JZ*, Lj = inf = inf{t > 0; R\ < L\}. (5.2) 

se[0,t] s€[0,i] 
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Finally, for k < I 6 ZU{— oo,+oo} and t € [0, oo], we consider the cr-field 



3k,l,t =<r{(oU), V s (j), N a (J), s G [0,t], I < j < k}. 



(5.3) 



Step 1. We first show that a.s., — r, < p 4 a.s., and that 770(1) = $i(Z l ), for 
some deterministic function <£>i, where 



Z l := (Co(i),iV s (j),K(i),i G {L> -l/2,...,i? pI + 1/2},sG [0,//]) . (5.4) 



The function <I> 4 can not easily be made explicit, see however Step 2 of the 
algorithm described in the Appendix [A"l 

It clearly suffices to treat the case i = 0. We consider the 60a; delimited by 
L\ and until they meet, i.e. until t = p°. Consider an avalanche process 
starting at some time T < ~p° with a given initial condition ip G E^ T . We wish 
to rebuild its value at time 0, thus the time goes now down on Figure 2 (see 
also Figure 4 below). 

Observe, having a look at Figure 2, that on the right and left of this box, the 
Bernoulli process £ f is vacant (see Lemma |4~5|) . so that due to our coupling, the 
avalanche process is also vacant, since it is always smaller (see Proposition ^. 4|) . 
Hence no interaction can go inside this box (from its left and right sides), since 
vacant sites cut the interaction: indeed, flocks falling outside this box can not 
make die flocks inside the box. 

Next, notice that the horizontal segments delimiting this box on the top side 
contain sites of the following type: 

(a) either vacant sites of the Bernoulli process, so that the avalanche process is 
also vacant at this site at this time (since it is always smaller); 

(b) either sites where the Bernoulli process becomes occupied because of a black 
mark, so that the avalanche process also becomes occupied at this time (because 
when the Bernoulli process is vacant, then the avalanche and Bernoulli processes 
both become occupied when they encounter a black mark) ; 

(c) either a grey mark (in the middle of an occupied zone of the Bernoulli 
process), so that at this site and at this time, the avalanche process is (or 
becomes) occupied. 

As a conclusion, the avalanche process is vacant on the left and right of the 
box, and the value of the avalanche process on the top horizontal segments of 
this box are determined, independently of its starting time T < —p° and initial 
condition cp G Eg T - We thus may rebuild the avalanche process {ijt' <f ')t£[p ,6\) 
and the obtained value ^'^(O) does not depend on T < —p° nor on ip G E^. 

We thus deduce that r > — p° and that 770(0) = 77^^(0). This value 770 (0) 
clearly depends only on the values of N, V, Ct in this box, and the fact that the 
Bernoulli process is vacant on the outside boundary of this box. We thus can 
say that 770(0) is a (deterministic) function of Z°. 

Step 2. We also observe that for i G Z, k < i < I, for T > 0, 



SU(k,l,T) := { 




(5.5) 
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Figure 4: Reconstruction of rj (Step 1). 




-5 -4 -3 -2 -1 1 2 3 4 5 



We explain here graphically how to obtain the values of r)o(i) for i £ {— 2, .., 1}. The avalanche 
process is represented on the right of each site, in grey. 

We start from the top of the box, that is at site 0. A flock appears here, so that becomes 

occupied (independently of its starting time T < — p° and initial condition ip £ ). 

Next, a grey mark appears at the site 3, which thus remains occupied (if it was already) or 

becomes occupied (if it was not), so that in any case, the site 3 is occupied at this time. 

At the same time, the sites 1 and 2 are vacant, since they are vacant for the Bernoulli process. 

Next, 1 and then —1 become occupied due to black marks. But then the flock at dies, and 

makes —1,0,1 become vacant. And so on... 

This way, we see that the site —1 is finally vacant at time 0, while the site is finally occupied. 
On the other hand, it is immediate that —2 and 1 are vacant, since the avalanche process is 
smaller than the Bernoulli process with our coupling. 

We could also see that r)o(2) = 1 and J7o(3) = here, but it is not possible to decide if 
?7o(— 3) = 1, because it could be killed by a flock dying at some site i < —6. 



and that Z 11 ^^.^^^)} is Cy/c^r-measurable. This is clear from Figure 2. 

Step 3. Next, we notice that point (a) (existence and uniqueness of an invariant 
distribution II for the avalanche process) follows immediately from Proposition 
13 - 1 1 provided we know that ti > — oo a.s. for all i £ Z. But we know from 
Step 1 that Ti > —p l . Lemma 1431 - fa) implies that for all i G Z, p l and p° are 
identically distributed. As a consequence, it suffices to show that p° < oo a.s., 
which follows from Proposition 14.41 (a) . 

Hence, II = C(j]o) is the unique invariant distribution of the avalanche process. 

Step 4. Point (c) (existence of a perfect simulation algorithm for (?7o(^))i=-i,...,i) 
is also immediate, see Appendix for an explicit simulation algorithm. Let I > 
be fixed. Due to Step 3, we know that rj ~ II, it thus suffices to simulate 
(perfectly) (r]o(j))je{~i,...,i}- This can be done by simulating first Z~ l ,...,Z l . 
This can be done due to Step 2, which says that Z l depends on Co, N, and V on 
an a.s. finite number of sites and on an a.s. finite time intervall. Next, it suffices 
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to compute T)q(i) — $>i(Z l ) for all i E {—/,..., I}, which can be done following 
the rules explained in Step 1, see also Figure 4. 

Step 5. We now check the mixing property anounced in point (d). Let thus 
n > 1 and k G Z be fixed. We consider the events (here [x] stands for the integer 
part of x G K) 

fi£ - < fe + [n/3] - 1/2} and ^ = > fc + [2n / 3 ] + 1/2} . (5.6) 

We deduce from Lemma 14.51(c) that on Q*, P^, < k + [n/3] - 1/2 for all 
i < k. We thus deduce from Step 2 that the family (Z 1 TLqi )i<k is £/-oo,fc+[n/3],oo~ 
measurable. Hence (jjo^ln 1 )i<fc is £_oo,fc+[n/3],°o- m easurable. By the same 
way, (r/o(i)ln2 )i>fc+n is <5fe+[2n/3], 00,00 -measurable. 

But of course, the two cr-fields £?_oo,fc+[n/3],oo an d £/fc+[2n/3],oo,oo are independent. 
Hence, for any A C {0, l}(-°°' fe ] , B C {0, l}l k + n >°°] , 

^ [(»ta(*))i<jfc g -4, (»?D(i))i>fc+n e P] 

-P [(vo(i))i<k G A] P [(Tto(i))i>fc+„ G S] 

<2P[(^) C ] + 2P[(^) C ], (5.7) 

so that 

^(.^fcju^,^ - n ( _ 00>fe] <g> n [fe+n)00) I < 2P((^) C ) + 2P((n£)°). (5.8) 

Using Lemma S3(a), we deduce that P[(fi r \) c ] = PfR^ > fc + [n/3] + 1/2] = 
PfpL > [n/3] + 1/2] and P[(fi» ) c ] = P[L^+" < fc + [2n/3] - 1/2] = P[i£, > 
n - [2n/3] + 1/2] < P[P^ > [n/3] + 1/2], since n - [2n/3] > [n/3]. Using finally 
Proposition 14.4K b) . obtain that for some constants 7 > 0, K > 0, 

pmm + mir] < 2P[p oo > [„/ 3 ] + 1/2] 

< 2e- 7 / 2 -T["/ 3 ]p[e^] < ife" 7 "/ 3 . (5.9) 
We thus obtain (L~3]) . setting q = e^/ 3 and C = 2if. 

Step 6. Finally, it remains to study the trend to equilibrium. For ip E E, 
t > 0, we denote by Tlf the law of the 95-avalanche process at time t. The main 
difficulty here is to obtain the trend to equilibrium for any initial datum ip. G E, 
since our coupling allows us a priori to deal only with initial data stochastically 
smaller than T. We thus have to introduce a final coupling, which mixes those of 
Lemma l2~3l and Proposition l3.ll We keep, however, all the notations introduced 
in this proof. We fix ip E E. 

Step 6.1. We consider a third IFPP W, independent of N, V, and Co- We 
consider, for T < using Lemma [2.31 the (£t, ip, N T , W T )-co\rp\ed Bernoulli 
processes (Ct,C?)te[T,o]- Recall that (C*)*g[t,o] ls the (<p, iV T )-Bernoulli process, 



27 



while (Cf T )ie[T,o] is the (<p, Mr)-Bernoulli process for some IFPP My. Notice 
that due to Lemma [2~3l (ii)- (iii) . we obtain that for any k > 0, any t € [T, 0], 

ifni(t,T,fc) := {<J(*) = C»(i); V<e [-*,*], a G[t,Q]} 

then P[Q. 1 (t,T,k)] > 1 - (2k + l) e - 2 ^- T \ (5.10) 

We finally consider the (ip,tf , Mt,V t )- coupled Bernoulli-avalanche processes 
(Ct,vI)te[T,o]- 

Step 6.2. We consider the event, for n > 0, t > 0, 



1/2, Li} > -I - n - 1/2, max < i 

ie{-z, •■•,*} 



(5.11) 



We know that on this event, riZ t,ip (i) = ?7o(i) for all i G {— I, as soon as 

tp G -E^. But we easily understand, using Step 6.1 and some arguments as in 
Step 1, that on f2i(— t, —2t, I + n + 1) PI Qo(n, t), we also have fjn 2t (i) = r)o(i) for 
all % G {-I, ...J}. 

Step 6.3. On the other hand, C(fjQ 2t ) = H% t and £(r) ) = II, so that we 
classically deduce that for all n > 0, 

| (n£) hM] - n hM] \ TV < 2P[(a (n, t)f u (n t (-t, -2t, i + n + i)) c ]. (5.12) 

Step 6.4. We obtain, using Proposition 14. 41 and Lemma l4~5l -(a). that 



P[(n (n,t)y] < p[3t,>i + n + 3/2 



P 



So' < -I - n - 3/2 



= -/ 



< 



2 > n + 3/2] + (2/ + l)P[p° > t] 



< 2e^ n E 



+ (21 + l)e~ 0t E 



< A(e^ n + (21 + l)e~ pt ) 



(5.13) 



for some constants A>0/3>0, 7>0. 

Step 6.5. Gathering ([STID]) . ([5T2"]) , fSTTg]) . we finally obtain that for any t > 0, 
any n > 1, 



|(nS)H,i]-n H)I] 



TV 



< 2A(e-~< n + (21 + l)e- f:it ) 
+2(2(1 + n+l) + l)e- 2t . 



(5.14) 



Choosing finally n = [t], we deduce that for a = min(l, 7, j3), we may find a 
constant K > such that for all t > 0, all I > 0, 



I (IF 



2t;[-M] 



n 



[-Ml I TV 



< K{l + l)e 



This yields (|1.3p . and concludes the proof. 



(5.15) 
□ 
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6 A related mean-field model 



This section, quite independent of the rest of the paper, is devoted to the brief 
study of a mean-field coagulation-fragmentation model related to the avalanche 
process. 

To obtain a process which preserves the total mass, we will slightly change our 
point of view: we assume that each edge of Z has a mass equal to 1. 
Consider a possible state 77 6 E of the avalanche process. Two neighbour edges, 
say (i — and + 1), are said to belong to the same particle if rj(i) = 1: 
the flock lying at i glues the two edges. For example, the edge (0, 1) belongs 
to a particle with mass 3 if r)(— 1) = rj(2) = and 77(0) = 77(1) = 1, or if 
77(0) = 77(3) = and 77(1) = 77(2) = 1. Similarly, (0, 1) belongs to a particle with 
mass 1 if and only if 77(0) = 77(1) = 0. 

Then we consider, for 77 e E and for k £ N, if it exists, 

. . number of particles with mass k in \—n, n] . . 

Ck (v) = lim 7T- n ' t 6 - 1 ) 

n— >oo Z77 + 1 

which represents the average number of particles with mass k per unit of length. 
Consider an avalanche process (r?t)t>o- Assume for a moment that for each t > 0, 
the successive masses of particles in 774 arc independent (which is intuitively far 
from being exact). Then, using the invariance by translation of the model, one 
would have, for k > 1, t > 0, 

Cfc(i) := Ck(rj t ) = \p [(0, 1) belongs to a particle with mass k in 774] (6.2) 
k 

The family (c(i)) t > = (cfe(i))t>o,fe>i would also satisfy J2k>i kc k(t) = 1 for all 
t > 0, and the following infinite system of differential equations: 

d Cl (t) = -2 Cl (t)+]T(fc-l)fc Cfe (t), (6.3) 



dt fc >i 



fe-i 



^-c k (t) = -2c k (t) - (k - l)c k (t) + ]— - y2 Ci (t)c k -i{t) (fc>2), 

dt m {c{t)) ^ 

where mo(c(t)) = X)fe>i c fe(0 i s the average number of particles per unit of 
length. Indeed, the first equation expresses that an isolated edge merges with 
its two neighbours with rate 1, while each time a flock falls on a particle with 
mass k, which happens at rate k — 1 (since a particle with mass k contains k 
edges and thus k — 1 sites) , an avalanche occurs and k new particles with mass 1 
appear. Next, the second equation expresses that particles with mass k become 
larger at rate 2 (when a flock falls on one of its two extremities), that particles 
with mass k disappear when they are subject to an avalanche (which happens at 
rate k — 1), and that particles with mass k do appear when a flock falls between 
a particle with mass i and a particle with mass k — i. This last event occurs at 
rate 1 , proportionnaly to the number (per unit of length) of pairs of neighbour 
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particles with masses i and k — i, which is exactly a(t)ck-i(t)/mo(c(t)). We use 
the abusive independence assumption when computing this last rate. 

We refer to Aldous [U Construction 5] for very similar considerations, without 
fragmentation, where the independence between neighbours really holds. 

The system (j6.3[) can be seen as a coagulation-fragmentation equation with 
constant coagulation rate K(i,j) = 2, with a splitting rate (from one particle 
with mass k into k particles with mass 1) F(k; 1, 1) = k — 1, the change of 
time l/mo(c(t)) lying in front of the coagulation term. Indeed, we could write, 
for example when k > 2, 

^c k (t) = -F(k;l,...,l)c k (t) (6.4) 

fc-i 

-c fc (t) K ( k > *) c i(<) + H k - i)ci(t)ck-i{t) . 

The term in brackets on the second line is the right-hand side member of the well- 
known Smoluchowski coagulation equation. See Aldous pQ , Laurengot-Mischler 
[14] for reviews on these types of equation. No result about trend to equilib- 
rium for such a model without detailed balance condition (here the coagulation 
is binary, which is not the case of fragmentation) seem to be available. See how- 
ever Fournier-Mischler [llj for partial results about a coagulation-fragmentation 
without balance condition in the case of binary fragmentation. 

We aim here to compute the steady state of this mean-field model, and to 
show numerically that it approximates closely the invariant distribution of the 
avalanche process. 

Proposition 6.1 The system of equations \6. 3\) admits a unique steady state 
c = (ck)k>i> in the sense that: c k > for all k > 1, ^2 k>1 kc k = lj and, with 
m o(c) = Efc>l c fc> 

2ci = y^(fc - i)fccfc, 

k>l 

- fc-1 

(k + l)c k = T-Vcift-j (fc>2). (6.5) 

m (c) f-{ 

This steady state is given by c k — a k g k ~ 1 2r k , where a\ = 1 and for k > 2, a k — 

YliZi a j a k-j, while g > is the unique solution ofJ2k>i a k(g/2) k = 1- 
We also have c\ = 1/2, X)fe>i k2c -k — 2, and rao(c) = 1/g. 

Proof Consider the sequence (a k ) k >i defined in the statement. Remark that 
for any x > 0, any g > 0, the sequence x\ = x, x k = ^Jl^ 1 XjX k -j is 
explicitely given by x k — a k x k g . 

Thus, (cfe)fc>i is a steady state of (|6.3p if and only if there exist x > and 5 > 
such that 



m (c(i)) 
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(i) V k > 1, c fe = a k x k g k - 1 , 

(ii) g = l/m (c), 

(iii) E fc >i fcc fc = 1> 

H^i(E t >i^-i). 

Points (i) and (ii) imply that necessarily, ^2 k>1 a,k(xg) k = 1. Thus g := xg is 
clearly uniquely defined, and satisfies < q < 1 (since ai = 1 and ci2 = 1/3 > 0). 
Next, using (iii), we deduce that g — 53fc>i ka kQ k is also uniquely defined (and 
finite, since q < 1 and since one easily checks recursively that ctfc < 1 for all 
fc > 1). Thus x = q/g is also uniquely defined. This shows that there exists 
at most one steady state. We next have to verify that these values for x and 
g imply point (iv). Using the definition of (a,k)k>i, we obtain that on the one 
hand, 

J2(k + l)a k q k = J2l k 12 a J ak -i = [E^ ) = *> ( 6 ' 6 ) 

k>2 k>2 j=l \j>l / 

while on the other hand, since a% = 1, 

J2( k + l Wq k = E kakqk + E akqk -2q = g + i-2q. (6.7) 

k>2 k>l k>l 

We obtain by this way g — 2q, so that x = 1/2, and thus c\ = 1/2. 

To conclude the proof, it suffices to check that J2k>i k 2 °k = 2 with the previous 

values for x and g. But again, we obtain on the one hand that 

Hk + l)a k q k = E fc2 °^ fe + E kakqk - 2< l = aJ2 k2 ° k +9-2Q, (6-8) 

fc>2 k>l k>l k>l 

while on the other hand, 

fc-l 

^2 K k + l)ak.q k = E fcgfe E a o ak -3 = E E k Q k ~ ja >k-j 

k>2 k>2 j=l j>l k>j + l 

= E a ^ J Ec? + l ^ ai = 2 ( E a ^ ( E lql ^ = 2 -9- ( 6 - 9 ) 

Hence gJ2k>i k2 °k + 9 ~ 2g = 2g, so that X)fc>i ^ 2c fc = 1 + 2<z/.g = 2. □ 

To conclude this section, let us give some numerical results. 
We obtain numerically, computing the values of a\,a2, ...,aiooooj and studying 
the function z i— > ^| 0000 ak(z/2) k , that g ~ 1.4458, with quite a good precision. 
We deduce then from Proposition 16. II that at equilibrium, the mean-field model 
(|6.3|) satisfies 

ci = 0.5, c 2 ~ 0.1204, c 3 ~ 0.04354, 
c 4 ~ 0.01679, c 5 ~ 0.006574, c 6 ~ 0.002582, 
Ek>iCk^ 0.6916, E fc >ifc 2 cfc = 2. 



31 



On the other hand, simulating 10 8 times the mass M of the particle containg 
the edge (0, 1) in the avalanche process 770 at equilibrium (see Appendix lA"]) we 
obtain the following Monte-Carlo approximations for Ck(rjo) '■= P[M = k}/k 

ci(ry ) ~ 0.499934, c 2 {m) - 0.12312, c 3 (?7 ) ~ 0.0422142, 
c 4 (r/ ) ~ 0.0161849, c 5 (r/ ) ^ 0.00648257, c 6 (r/ ) ~ 0.00263739, 
E fe >i c fc (r?o) - 0.692419, J2 k >i k 2 c k ( V o) - 1.99979. 

It appears clearly that the two sets of values are very similar, even if numerical 
computations indicate that no equality holds, except maybe concerning C\ and 
y~] k k 2 c k . We have no explanation for this phenomenon. It might indicate that 
correlations between the masses of successive clusters are nearly insignificant. 
We have no proof that the mean-field model is the (very fast) limit, in some 
asymptotic regime, of the avalanche process. 

A Appendix 

Let us now write down the simulation algorithm, which we deduce from Sec- 
tions [3] and 2J For I > 0, the algorithm below simulates a random variable 
(Vo(i))ie[-i,i] with law IIr_; ) ;i, where II is the unique invariant distribution of 
the avalanche process. The idea is to simulate N, £ and V in an a.s. finite 
random space-time domain and then to reconstruct r/ following the graphical 
construction 12.51 

We construct a random process ( n (k) containing the values of £ at some random 
times T n (times of jumps of N and V in a finite growing spatial domain [£ n , r n ]) 
and an additionnal information: roughly, 
Cn(fe)=0ifCf„(fc) = 0; 

(n(k) = 2 if (k) = 1 and (k, T n ) belongs to the box delimited by the contour 
processes; 

(n(k) = 1 if Cf (k) = 1 and (k,T n ) is outside the box. 

We invite the reader to have a look to Figures 2 (for Step 1) and 4 (for Step 2) 
while reading the simulation algorithm below. We will say the box for the box 
delimited by the contour processes. 

Simulation Algorithm for Or_i ; n 

Step 0: Initialization. 

Simulate the initial Bernoulli configuration Co(^) for k £ [?o,ro], where Iq (resp. 
ro) is the first vacant site on the left (resp. right) of — I (resp. I). 
If (o(k) = for all k 6 [— /, /], set rjo(k) = for all k € [— /, I], and stop here. 
Else, set Co(^) = 2Co(^) for k G [lo, r ], set n — 0, and proceed to Step 1. 

Initially, all the sites in [Iq + 1, ro — 1] are in the box. We thus assign the value 
to vacant sites, and the value 2 to occupied sites. 

Step 1: Simulation of black/grey marks and contour processes 

Set n = n + 1 . 
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Choose i n uniformly in [£ n -i,r n -i] representing the involved site. 

Choose m n ~ Ber(l/2), here m n — for a black mark, m n — 1 for a grey mark. 

• If m n = 0, and ( n -i(i n ) > 1 then we set ( n (i n ) = 0. 

// the site i n is occupied, it becomes vacant due to a black mark. 

• If m n = 1 and ( n -i(i n ) = 2; 

set Cn(i n ) = 1 if C„-i(i„ - 1) = Cn-l(*n + 1) = and if Vfc G [^„_i,i„ - 
1], Cn-i(k) < 1 or Vfc G [i n + l,r„_i], C„-i(fc) < 1; 

The site i n remains occupied but leaves the box due to a grey mark, because 
its two neighbors are vacant, and because it is on the boundary of the box. 

otherwise, set ( n (in) = 2. 

The site i n remains occupied and in the box, because either the site i n is 
in the strict interior of the box or one of its neighbors is occupied. 

• If m n = and ( n -i(in) — then we consider the connected component 
I n of occupied sites (plus i n ) around i n (at time n — 1). 

If Vfc G InU[i„ + l,r„-i], Cn-i(fc) < 1 orifVfc G 7„U[£„_i,i„-l], ( n {k) < 
1, then, for all k e I n , set ( n (i n ) = 1- 

Otherwise, set ( n +i{k) — 2 for all k G I n . 

The site i n becomes occupied due to a black mark. Then its whole connected 
component of occupied sites joins the box, except if all these sites were 
outside the box at time n — 1 . 

• Set ( n (k) = Cn-i(fc) for all sites l n -\ < k < r„_i of which the value (at 
time n) has not been defined yet. 

We update all other sites. Observe that we have not considered the case 
m n = G {0,1}: grey marks have no effect on vacant sites in 

the Bernoulli process, and cannot increase the number of sites in the box. 

If Cn( r n-i) < 1) se t T n — r n -i- Else, consider s r n ~ Geo(l/2) (i.e. P[s r n = k] = 
(l/2) fe for k > 1), set r n = r n -\ + s n , ( n (k) = 2 for k G [r„_i,r„ - 1], and 
Cn{r n ) = 0. 

We extend the box if necessary, i.e. when Cn( r n-i) is occupied and in the box. 
We thus extend the Bernoulli process to the right until we meet a vacant site, 
at some site r n . Then the (occupied) sites k G [r„_i + 1, r n — 1] are in the box. 

If Cn(^n-i) < 1) set i n = £ n -i- Else, consider s n ~ Geo(l/2), set £ n = £ n -i — s n , 

(n(k) = 2 for fc G [l n + ljn-l], and Cn(^n) = 0. 

Here we use the same arguments on the left of the domain. 

Check {l n < k < r„,(^„(fc) = 2}. If this set is non empty then repeat Step 1. 
Otherwise, set T — n and proceed to Step 2. 
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If all the sites have a value equal to or 1, this means that the contour processes 
have met. 

Step 2: Deduction of the avalanche invariant realization. 

Start with fj T (k) = for all l T < k < r T . Then for all 1 < n < T, define 
recursively (for n decreasing from T to 1) r} n -i(k), for all £ n -i <k< r n _i, by 

• Vn-i(k) = if m„ = and if k belongs to the connected component of 
occupied sites of i n (in rj n ). 

Black marks kill connected components of occupied sites. 

• rjn-i(k) = 1 if k = i n , if m n = 1 and if ( n -i(i n ) > 1. 

Grey marks make appear flocks when the Bernoulli process is occupied. 

• Vn-i(k) = 1 if k = i n , if m n = and if f} n (i n ) = C n (i n ) = 0. 
Black marks make appear flocks at vacant sites. 

• Vn-i(k) = rj n {k) for all sites k e [£ n -i,r n -i\ for which the value r/ n -i(k) 
has not been defined yet. 

Conclusion. 

Then {rj (k),k e [— 1,1]} is distributed according to U[-i^. 

Remark that the number of iterations of Step 1 is finite due to our results: when 
the contour processes L~ l and R l meet, there are no more sites with value 2. 

Alternative 

We finally propose another version of Step 1. The advantage is that the number 
of involved sites is much smaller. The idea is to take better advantage of the 
so-called grey marks : we will keep track only of what may really be needed to 
reconstruct {770(^)7^ £ [— M]}. 
Step 1'. 

Set n = n + 1. Choose i n uniformly in [£ n -i, r n -i]- Choose m„ <~ Ber(l/2), 

• If m n = 0, and ( n -i(i n ) > 1 then we set ( n {i n ) = 0. 

• If m n = 1 and Cn-i(^n) = 2; we set ( n {in) = 1 as soon as ( n -i(i n — 1) < 1 
or Cn-i(«™ + 1) < 1; else we let Q n {i n ) = 2. 

• If m n = and Cn-i(*n) = then we consider 

/„+ = {ke [i n + l,r„_i],Cn-i(fc) = 2, Vl„ < j < fc,Cn-l(j) > 1}, 

In- = {kE [e n -l,in - l],Cn-l(fe) = 2, Vfc < j < t„,Cn-l(j) > !}• 

If /„+ = /„- = 0, we set Cn(in) = 1- 

Else, we set ( n (k) = 2 for all k G [min(7 n — ), max(7 n +)]. 
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• Set Cn(k) = Cn-i(fc) for all sites £ n -i < k < r„_i of which the value (at 
time n) has not been defined yet. 

If there is k £ [£ n -i> i"n-t] such that Cn(k) — 2 and Cn(i) > 1 f° r au < J < 
k, consider s l n ~ Geo(l/2), set £„ = t n _ x -sl n , ( n (j) = 1 for j 6 [4 + 1, 4_i - 1] 

and Cn(4) = 0. 
Otherwise, set l n = £„_i. 

Act symmetrically on the right of the domain. 

Check {£ n < k < r n ,C, n (k) = 2}. If this set is non empty then repeat Step 1'. 
Otherwise, set T = n and proceed to Step 2. 

Let us emphasize the differences between Step 1 and Step 1'. The domain is 
extended only if there is k with Cn(k) = 2 in the connected component touching 
the boundary. The first and fourth switching rules of ( n (k) are left unchanged. 
In the second one, it is made easier to set Cn(k) from 2 to 1. In the third one, 
a much smaller set of indices is switched from 1 to 2. 

One may understand that changing Step 1 into Step 1' does not change the law 
of the final values {vo(k))ke[-l,l]- It defines, in some sense, a much more precise 
contour process (leftmost and rightmost sites k with ( n {k) = 2) than the one 
defined in Section [4j we numerically observe that using Step 1', the algorithm 
is 15 times faster than when using Step 1. We hope to take advantage of this 
idea to generalize our methods to more general particle systems. 
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